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o: 

O; 

■ Abstract 

|J^ \ Linear response theories in the continuum capable of describing continuum spectra 

Qv ' and dynamical correlations are presented. Our formulation is essentially the same as the 
continuum random-phase approximation (RPA) but suitable for uniform grid representa- 

^7j ! tion in the three-dimensional (3D) Cartesian coordinate assuming no spatial symmetry. 

^f^ ' Effects of the continuum are taken into account either by solving equations iteratively 

CN ■ with a retarded Green's function or by an absorbing boundary condition. The methods 

^^ \ are applied to giant resonances in a deformed nucleus ^^C. 

o: 

O ; 1 Introduction 

"*7^ ! Recent developments in radioactive beam techniques enable us to access nuclei near and 

Q ' beyond drip-lines. The drip-line nuclei are weakly bound finite fermion systems. One 

^ . would naturally expect that the continuum should be taken into account explicitly in 

,. \ description of their structure and reaction. 
^ ' Bound solutions of the Schrodinger equation for the Hamiltonian H = T + V 

'k/ \ 

\:\- {E-T)^ = V^, (1) 

c^ \ 

can be formally written as 

^ = T^^y^- (2) 

E-T ^ ' 

Here, a Green's function {E — T)^^ can be uniquely defined since the operator E — T 
is a negative-definite operator. Bound states are characterized by discrete spectra, En 
{n = 0,1,2,- ■■). 

On the other hand, for continuum states, since the E — T has zero eigenvalues, we 
need to modify Eq. (|2|) into the Lippmann-Schwinger equation. First, we need to add 
a zero-eigenvalue solution of the E — T, ^q. Furthermore, in order to uniquely define 
the Green's function {E — T)~^, it is necessary to specify a boundary condition. For 
usual physical situations, it is natural to adopt an outgoing boundary condition (OBC) 
for scattering waves, 

h — 1 + irj 
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For two-body systems, one may easily separate relative coordinates from those of the 
center of mass. Then, the OBC becomes rather trivial because this is essentially a 
problem of a single degree of freedom. For many-body systems, however, it can be very 
complicated to settle the OBC for many degrees of freedom. We focus our discussion 
upon treatment of this problem. 

The Skyrme Hartree-Fock (HF) theory in the three-dimensional (3D) Cartesian coor- 
dinate representation Q has been extensively applied to study of ground-state properties 
of relatively heavy nuclei. To investigate their excited states, a straightforward exten- 
sion is either the time-dependent Hartree-Fock (TDHF) or the linear response calculation 
based on the HF ground state. For spherical nuclei, since the HF problem has a single 
degree of freedom (the radial coordinate from the center of mass), one can easily treat 
the boundary condition. Therefore, the inclusion of the single-particle continuum for 
particle-hole (p-h) excitations has been achieved by using a method, so called, contin- 
uum random-phase approximation (RPA)|^]. The continuum RPA combined with the 
Skyrme Hartree-Fock (HF) theory has been extensively utilized to study giant resonances 
(GR) in spherical nuclei||3|, Q|. However, it is not trivial to extend the method to deformed 
nuclei. 

Since the Hamiltonian in the Skyrme HF theory is almost diagonal in coordinate 
representation, a grid representation in the coordinate space provides an economical 
description. The main issue is then how to incorporate the boundary condition in the 
3D uniform grid representation. We shall present prescriptions to treat the full three 
dimensional (3D) continuum in the RPA level. 

Recently, we have investigated photoabsorption cross sections of small molecules with 
essentially the same techniques |^]. In this paper, we report the applications to nuclear 
giant resonances with a simplified Skyrme Hamiltonian. 

2 Linear response in the continuum 

2.1 Outgoing boundary condition (OBC) in the 3D space 

In the linear response theory based on the Skyrme HF, exact treatment of the single- 
particle continuum is possible utilizing a single-particle Green's function. For spherical 
systems, the Green's function can easily be constructed by making a multipole expansion 
and discretizing the radial coordinate. This is an essential part of the continuum RPA 
method. In this section we present a method to construct a single-particle Green's 
function in the 3D grid representation for a system without any spatial symmetry. 

The linear response theory is formulated most conveniently using a retarded density- 
density correlation function Q. 

"'"'-' ^ /''-"'-""('■'^''■»). (^) 

in(r,t;r',t') = 9(t - t')(0|[p(r,t),p(r',t')l|0)' P) 

In the RPA of retaining only ring diagrams, H(r, r'; uj) is constructed from independent- 
particle density-density correlation function, Hq, which is defined by identifying the state 
|0) in Eq. (g) with the unperturbed HF ground state and assuming that the density 
operator is evolved in time with the static HF Hamiltonian. Imn(a;) is proportional to 
the excitation probability from the ground state to a state of excitation energy uj. 
The Hq is written in a formQ 

A 

Ho(r,r';^) = 5^<A.(r){(GW(r,r';(e,-^)*))* + G(+)(r,r';e, + .;)}<A.(r'), (6) 



assuming that the occupied states have real wave functions (pi and eigenenergies e^. The 
G'^'^' is a single-particle Green's function with OBC, defined by 



G^+Hr,r';E) = {r\(E + ^V^-VnF{r)+ir]) \r'), 
\ 2m J 



(7) 



where Vhf is the HF potential which is not necessarily spherical. It is difficult to explicitly 
construct the G^^'{r,r'; E) in the 3D space. For bound states {E < 0), instead, we can 
calculate actions of G^~^^ (E) on a state $ by solving a differential equation 

E + ^V2 - yHFCr)) ^(r) = $(r), (8) 

where ^ = G^^'{E)^. However, as we have mentioned in Sec. |l|, we have to specify the 
OBC for the continuum states with E > 0. 

For this task, we decompose the HF potential into a spherical part Vq and a short- 
range deformed part Vo (Vhf = Vq + Vr,). We first construct a Green's function, 
g'1)~^\e) = {E - Ho + i-ny^, for the spherical Hamiltonian, Hq = -l/2mV'^ + Vq, 
by using the usual multipole expansion technique Q. Next, we solve an equation 

il-Gi;\E)Vj,)^ = Gi;\E)<^, (9) 

to calculate * = G'^'^^E)^ = {E - Hq - Vd + iry)~^$. Eq. (|) can be derived from an 
identity for the Green's function 

G(+)(E) = Gi~'\E) + Gi-'\E)VBG^+\E). (10) 

In this way, we fix an outgoing asymptotic behavior of the Green's function G^^'{E). 
The details of numerical procedure is described in our recent paper[^. 

2.2 Absorbing boundary condition (ABC) equivalent to 
OBC 

In this section, we present another method to simulate the OBC, which we call "absorbing 
boundary condition" (ABC)|^, ^. The Green's function with OBC is written as Eq. (^) 
where +ir] is an infinitesimal imaginary quantity (r/ > 0). Now, we allow this imaginary 
part to depend on coordinate and to be finite, +iW{r). W{r) is taken to be positive far 
outside the system (at large r) and zero elsewhere. This means that, for E > 0, the wave 
number k has a positive imaginary part +i7(r) at large r. The outgoing wave defines 
its asymptotic behavior as 

V'(+) ~ /(J^)^, atr^oo. (11) 

r 

Thus, with the complex potential —iW{r), the outgoing wave is going to damp as 

pikr—'yr 

V-^^^ ~ f{n) , at r ^ oo, (12) 

while the incoming wave is going to diverge. Therefore, if we impose the vanishing 
boundary condition at large r, only outgoing waves are allowed. Since the complex 
potential takes care of the boundary condition, we can directly solve Eq. (|8|) with 
^ I boundary = 0, instead of Constructing the OBC by solving Eq. (^. In this sense, the 



ABC is easier to handle than the OBC in Sec 2.1. 



The remaining task is how to construct a good absorbing boundary potential. When 
an outgoing wave hits the complex potential —iW{r), a part of the wave is absorbed 
and another part is reflected back. Since the reflected part causes signiflcant inaccura- 
cies, the complex boundary potential must be strong enough to absorb the whole wave 
and simultaneously gentle enough to minimize the reflection. We adopt an absorbing 
potential of linear dependence on the radial coordinate, 

w/ \ /O' forr <i?, 

W{r) = ^W^r_^^ for R<r<R + At. ^^^^ 

In order to minimize the spurious reflection, the height Wq and width Ar should satisfy 
a condition|^, |5| 

20 ^= < \Wo\ < —ArVs^E^/^. (14) 

ArySm 10 

An advantage of the ABC is simplicity of its numerical calculation. We may solve 
scattering problems in the same way as we do for bound states. Another advantage is 
that the ABC easily implements a real-time calculation |0, |5| . Instead of taking an energy 
representation, we may calculate time evolution of the TDHF states directly. We have 
demonstrated that the real-time method with the ABC can properly take account of 
the continuum in the linear response calculations |^]. In real-time calculations, it is very 
difficult to handle the OBC explicitly. The disadvantage is that the ABC cannot handle 
low-energy escaping particles and a long-range Coulomb potential properly. Therefore, 
the outgoing protons feel a finite-range Coulomb potential, Ze/r only at r < i?. 

3 Giant resonances in the continuum 
3.1 Giant monopole resonances 

We apply the methods to isoscalar giant resonances (GMR) in ^^C. We use an energy 
functional given by the "BKN Hamiltonian" ||lO| ] : 

r A 
^^P^ = h^ ^El^'^*Wl' + i*o^W' + l^*3P(r) 

^ J J |r-r'| 

+ ^JdrJdr'pir)j-^^pir'), (15) 

where spin-isospin degeneracy (each nucleon with a charge e/2) is assumed, p{r) = 

The calculated HF ground state of ^^C has an oblate shape. Ratio of the minor 
and major axises is about two to three. We carry out three kinds of linear response 
calculations with different boundary conditions; (1) Box boundary condition (BBC) in a 
3D coordinate space of radius R = 18 fm, (2) Absorbing boundary condition (ABC) in 
a space of i? = 8 fm plus Ar = 10 fm, and (3) Outgoing boundary condition (OBC) in 
a space of i? = 8 fm. The BBC simply means a vanishing boundary condition, namely 
ip{r) = at r > R fm. A mesh spacing is taken as Ax = Ay = Az = 1 fm. In order to 
produce a smooth response curve with BBC, we adopt a complex energy E + ir/2 with 
r = 1 MeV. First, we apply an external field of Mqq = r^ to see a monopole response. 



Results of the calculation are shown in Fig. 3.1. The BBC provides a wrong result even 
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(a) BBC with r=1 MeV 
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Figure 1: Monopole strength distribution, 7rp{E)\{E\Moo\0)\'^, for ^^C. 
The sohd hues are the RPA strengths while the dashed hues indicate the 
unperturbed ones. The same smoothing parameter F = 1 MeV is used for 
three different calculations, (a) BBC (Box boundary condition), (b) ABC, 
and (c) OBC. See text for details. 



though we have used a large box of i? = 18 fni. The result of ABC is very similar to 
that of OBC, except for small discrepancies in low-energy part (E < 19 MeV). Since the 
highest occupied single-particle energy is —13.9 MeV, energies of outgoing particles are 
less than 6 MeV in this energy region. Actually, using Ar = 10 fm, it is impossible to 
satisfy the condition, Eq. (|T^, for these low-energy outgoing particles. We need a larger 
model space (larger Ar) to satisfy the condition. Eq. ( |l4|) is satisfied at higher energy, 
E — 13.9 MeV > 10 MeV. The ABC and OBC calculations seem to indicate two peaks 
at -E Rs 20 MeV (FWHMrs 5 MeV) and at E ^ m MeV (FWHM~ 15 MeV). 

Next, let us discuss isoscalar giant quadrupole resonances (GQR). The results are 
shown in Fig. 3.1. For GQR, deformation splitting is well-known. The K = 2 quadrupole 
strength is well localized around 16 MeV, for which the BBC calculation well reproduces 
the OBC result. On the other hand, K = {) quadrupole strength is located at higher 
energy and possesses a wider width. The single peak at i? ~ 21 MeV splits into three 
peaks in the BBC calculation, although the average energy over the three peaks be- 
comes a correct value. Fig. ^T] (b) also represents the monopole strength (dotted line). 
The results seem to suggest a deformation mixing between the monopole and K = Q 
quadrupole excitations around £" ~ 20 MeV. 
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Figure 2: Quadrupole strength distribution for ^^C. The sohd hues are the 
K = quadrupole strength, 7rp{E)\{E\Q2o\0)\'^, while the dashed lines indi- 
cate the K = 2, TTp{E)\{E\Q22\0)\'^. A panel (a) is the result of calculation 
with BBC, and (b) is the one with OBC The ABC gives the result almost 
identical to (b). The smoothing parameter F = 1 MeV is used. The dotted 
line in a panel (b) indicates a monopole strength distribution in arbitrary 
units. 



4 Conclusions 

Outgoing boundary condition (OBC) is one of the most important ingredients to treat 
the continuum. In order to settle the OBC in the calculation on the 3D mesh space, 
we have presented two methods. The exact treatment can be done by using a Green's 
function technique, making resort to an iterative numerical procedure to solve Eq. (^). 
Another method is to use an absorbing boundary potential. Carefully choosing the 
complex potential, we may simulate the OBC with the same procedure as we do for 
bound states. We show some test applications to isoscalar giant resonances with the 
BKN Hamiltonian. The conventional box boundary condition may give wrong results for 
the continuum response even if adopting a large box (R ~ 20 fm). The two methods are 
complementary in a sense that the first method treats the OBC exactly and the second 
becomes a very efficient method for high-energy outgoing particles. We are preparing 
applications of present methods to response and reactions in drip-line nuclei. 

References 

[1] H. Flocard, S. E. Koonin, and M. S. Weiss, Phys. Rev. C 17, 1682 (1978); 
K. T. R. Davies, H. Flocard, S. Krieger, and M. S. Weiss, Nud. Phys. A 342, 
111 (1980). 

[2] S. Shlomo and G. Bertsch, Nud. Phys. A A243, 507 (1975). 

[3] N. van Giai and H. Sagawa, Nud. Phys. A A371, 1 (1981). 



6 



[4] I. Hamamoto, H. Sagawa, and X. Z. Zhang, Phys. Rev. C 53, 765 (1996); 
I. Hamamoto and H. Sagawa, Phys. Rev. C 53, R1492 (1996); Phys. Rev. C 54, 
2369 (1996). 

[5] T. Nakatsukasa and K. Yabana, J. Chem. Phys. 114, 2550 (2001). 

[6] A. L. Fetter and J. D. Walecka, Quantum Theory of Many -Body Systems (McGraw- 
Hill, New York 1971). 

[7] I. Hamamoto and B. R. Mottelson, J. Phys. Soc. Japan Suppl. 44, 368 (1978). 

[8] R. Kosloff and D. KoslofF, J. Comput. Phys. 63, 363 (1986). 

[9] M. S. Child, Mol. Phys. 72, 89 (1991). 

[10] P. Bonche, S. Koonin, and J. W. Negele, Phys. Rev. C 13, 1226 (1976). 



